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ABSTRACT 

Aims. We use the deepest existing mid- and far-infrared observations (reaching ~3 mjy at 70yum) obtained with Spitzer in the 
Great Observatories Origins Deep Survey (GOODS) and Far Infrared Deep Extragalactic Legacy survey (FIDEL) fields to derive 
the evolution of the rest-frame 15 yum, 35 /urn, and total infrared luminosity functions of galaxies spanning z < 1.3. We thereby 
quantify the fractional contribution of infrared luminous galaxies to the comoving star formation rate density over this redshift range. 
In comparison with previous studies, the present one takes advantage of deep 70 pm observations that provide a more robust infrared 
luminosity indicator than 24 yum affected by the emission of PAHs at high redshift (z ~ 1), and we use several independent fields to 
control cosmic variance. 

Methods. We used a new extraction technique based on the well-determined positions of galaxies at shorter wavelengths to extract 
the 24 and 70 pm flux densities of galaxies. It is found that sources separated by a minimum of 0.5xFWHM are deblended by this 
technique, which facilitates multi-wavelength associations of counterparts. Using a combination of photometric and spectroscopic 
redshifts that exist for ~80% of the sources in our sample, we are able to estimate the rest-frame luminosities of galaxies at 15yum 
and 35 pm. By complementing direct detections with a careful stacking analysis, we measured the mid- and far-infrared luminosity 
functions of galaxies over a factor ~100 in luminosity (10" L G <~ Lir <~ 10 13 L Q ) at z < 1.3. A stacking analysis was performed 
to validate the bolometric corrections and to compute comoving star-formation rate densities in three redshift bins 0.4 < z < 0.7, 
0.7 <z< 1.0 and, 1.0 < z < 1.3. 

Results. We find that the average infrared spectral energy distribution of galaxies over the last 2/3 of the cosmic time is consistent 
with that of local galaxies, although individual sources do present significant scatter. We also measured both the bright and faint ends 
of the infrared luminosity functions and find no evidence for a change in the slope of the double power law used to characterize the 
luminosity function. The redshift evolution of infrared luminous galaxies is consistent with pure luminosity evolution proportional to 
(1 +z) 3 ' 6±04 up to z ~ 1.3. We do not find evidence of differential evolution between LIRGs and ULIRGs up to z ~ 1.3, in contrast with 
previous claims. The comoving number density of infrared luminous galaxies has increased by a factor of ~ 100 between < z < 1. 
By z ~ 1.0, LIRGs produce half of the total comoving infrared luminosity density. 
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1. Introduction 

Constraining the star formation history of galaxies as a function 
of redshift is a key to understanding galaxy formation. To 
estimate the ongoing star formation rate (SFR) of galaxies one 
needs to measure direct and re-radiated stellar emission of 
young stars which correspond to the rest-frame UV and IR light, 
respectively. Due to the difficulties and the only recent possi- 
bility to obtain deep infrared observations, first studies which 
have estimated SFR evolution with redshift ha ve only used 
UV o bservations corrected for dust extinction (Mada uet al.l 
I1999I) . However dust extinction correction suffers from various 
uncertainties and it is well known that local galaxies harboring 
strong dusty star formation are opaque to UV radiation (Buat 
et al. 2005). Hence to obtain a complete understanding of these 
dusty star forming galaxies, deep infrared observations are 
required. 

Taking advantage of infrared capabilities such as the 



Infrared Space Observatory (ISO), several teams have confirmed 
the import ance of infrared galaxies for the understanding of SFR 
evolu ti on fElbaz et alj|1999t lAussel et alj|1999t [Chary & Elbaz] 
I 2001 I : iFranceschini et all I200U IXuetalJ 120011: lElbazetalj 
l2002l:lMetcalfe et alJl2003l:lLagache et alj|2004l) . Indeed they all 
found that the contribution of luminous infrared galaxies to the 
SFR density increases with redshift up to z ~ 1 and that at high 
redshift the bulk of the SFR density occurs in luminous infrared 
galaxies (i.e 10 U L G <LIRGs< 10 12 L G , 10 12 L G <ULIRGs) 
which are relatively rare in the local Universe. 

The Spitzer Space Telescope made it possible to obtain 
a complete census of the SFR density evolution up to z ~ 1 
using deep 24 u rn Multiband Imaging Phot ometer (M IPS; 
Riek e et all 2004 ) observations. Many authors dL e Floc'h et all 
l2005UCharv & Elbazll200lHLagache et alJl^llXtiEoOOj) have 
studied the evolution of the rest-frame 15 pm or 8 pm luminosity 
function (LF) ba sed on MIPS 24 pm dat a. Then using empirical 
SED libraries dCharv & Elbazl 1200 U lLagache et all 120031: 
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iDale & Heloul 120021) to derive the total infrared luminosities 
of galaxies from these monochromatic luminosities, they were 
able to constrain the evolution of the total infrared LF. They 
all found a strong evolution with cosmic time of the relative 
contributions of normal, LIRG and ULIRG galaxies to the SFR 
density. At z ~ 0, the SFR density is dominated by quiescent 
galaxies whereas at z ~ 1 it is dominated by LIRGs. They 
also found that the contribution of ULIRGs to the SFR density 
rises steeply from z — to z = 1 but still remains negligible at 
z ~ 1 (i.e. less than 10% ). More recently, studies using MIPS 24 
pm dCaputi et al.ll2007h data have extended previous works to 
higher redshift and found that ULIRGs may dominate the SFR 
density at z ~ 2. 

All these studies, based on 24 pm observations, are 
strongly dependent o n the SEP library (ICharv &Elba3 l200ll 
lLagache et ail 120031; IDale & Heloul 120021) used to derive the 
total infrared luminosity of galaxies. Moreover significant 
uncertainties on the bolometric correction are implied by the 
redshifting of the strong PAH and silicate emission into the 24 
pm bandwidth. 

At z ~ 0, ICalzetti et ail d2007l) have proven that the rest- 
frame 24 pm luminosity is a rob ust SFR estima t or. Mo reover 
using deep radio observations Ap pleton et al.l (120041) have 
proven that observed 70 pm flux was a more robust SFR 
estimator than the observed 24 pm flux. Hence using the 
observed MIPS 70 pm data for galaxies below z ~ 1.3, as in 
the present study, will strongly reduce the bolometric correction 
uncertainties since this emission corresponds to the rest-frame 
35 pm luminosity at z ~ 1 which is produced by dust of about 
the same temperature (T ~ 40 K) as that responsible for the 24 
pm emission in local galaxies. The aim of the present paper is 
to take advantage of the deepest 70 pm observations currently 
available to strongly constrain the bright end slope of the 35 pm 
and the total infrared LFs up to z ~ 1.3. These far-infrared data, 
as part of the Far Infrared Deep Extragalactic Survey (FIDEL) 
and the Great Observatories Origin Deep Survey (GOODS), 
cover 4 fields (Extended Groth Strip, Extended Chandra Deep 
Field South, GOODS-N and, GOODS-S) which are wide 
enough (~ 1400 total arcmin 2 ) and deep enough (~ 3mJy at 
70 /im) to obtain a good constraint on the bright end slope of 
these LFs. By complementing direct 70 /im detections with a 
careful stacking analysis, we also constrain the faint end slope 
of these LFs. Finally, by studying the redshift evolution of the 
total infrared LF, we derive the evolution of the SFR density and 
the relative contributions of normal, LIRG and ULIRG galaxies 
to the total comoving infrared luminosity density as a function 
of redshift. 

The layout of the paper is as follows: In Section [2] 
we present the mid- and far-infrared observations as well as 
all the ancillary data used in this study. We also present in this 
section the new extraction technique used to extract infrared 
sources and the simulations made to characterize the limits of 
our infrared catalogs. In Section[3] we present the rest-frame 15 
pm LF inferred in three redshift bins center at z — 0.55, 0.85 
and, 1.15 using MIPS 24 /mi data. The rest-frame 35 pm LF 
inferred in these three redshift bins using the MIPS 70 /im data 
is presented in the Section [4] In this Section we also measure, 
using a stacking analysis, the observed 24/70 pm correlation 
which enables us to constrain the faint end slope of the rest 
frame 35 pm LF. In Section [5] we present the total infrared LF 
function and its evolution up to z ~ 1.3. Finally in Section 
[6] we discuss the evolution of the rest-frame 35 fim and total 
infrared LF in terms of the cosmic SFR history and the relative 



contributions of normal and luminous infrared galaxies. 



Ho 



Throughout this paper we will use a cosmology with 
= 71 km s' 1 Mpc-\Q. A = 0.73, Q M = 0.27. 



2. The Data 

2.1. Infrared Imaging 

2.1.1. EGS 

The Extended Groth Strip field (14*17 m , +52°30') was ob- 
served at 24 and 70 fim with the MIPS instrument on board 
the Spitzer Space Telescope as part of the Far Infrared Deep 
Extragalactic Legacy program (FIDEL, PI: Dickinson). These 
observations cover a total area of 900 arcmin 2 using the scan 
mode AOT of Spitzer (see Spitzer observer's manual for more 
information). 

The data were coadded using the MOPEX mosaicing soft- 
ware provided by the Spitzer Science Center. The final 24 fim 
mosaic has a pixel scale of 1.2"/pixel and an effective integra- 
tion time per sky pixel of ~ 14000s. The final 70 fim mosaic 
has a pixel scale of 4.0" /pixel and an effective integration time 
per sky pixel of ~ 7 200s. For both wavelengths the calibra- 
tion factor used to convert digital units (DN/s) to flux (MJy/sr) 
was taken from the MIPS data handbook (see SSC MIPS Data 
Handbook v3.3.1 section 3.7.2). At these wavelengths the point 
spread functions (PSF) are characterized by a FWHM of ~ 5.9" 
and ~ 18" at 24 and 70 fim observations, respectively. These 
FWHM were measured on the empirical PSF used to perform 
our PSF fitting analysis (see Section l2~2l >. 

The 24 and 70 fim observations reach a point source sen- 
sitivity of ~ 50 /iJy and ~ 3 000 /iJy respectively. These point 
source sensitivities include effects of background and confusion 
noise and were derived using extensive simulations (see Section 
12. 21 for a complete definition). For both wavelengths, source ex- 
traction was first performed with a PSF fitting technique using 
the positions of IRAC sources as priors (see Section [2T2l . This 
procedure leads to a total of 8457 sources at 24 fim (> 50/Jy) 
and 669 sources at 70 fim (> 3 mJy). 

2.1.2. GOODS-North & South 

Infrared imaging of the GOODS-N field (12 A 36 m , +62°14') 
was obtained at 24 pm as part of the GOODS Legacy program 
(PI: Dickinson ). The 70 pm data were obtained as part of the 
GO-3325 (PI: iFraver et al.l l2006bl) program for the 10' x 10' 
central region and as part of the FIDEL legacy program for 
the northeast and southwest regions in order to cover the 
entire GOODS-N area (i.e 15' x 10'). Infrared imaging of the 
GOODS-S field (3 h 32 m , -27°48') was obtained at 24 and 70 
pm as part of the GOODS legacy program (PI: Dickinson) and 
the GO-20147 (PI: Frayer) program respectively. We note that 
even if the GOODS-S 24 pm field covers a 10' x 18' area, the 
deep 70 pm observations only cover a 10' x 10' area. Thus the 
GOODS-N and the GOODS-S fields cover respectively a total 
area of 194 arcmin 2 and 91 arcmin 2 . 

Final mosaics, obtained using MOPEX, have the same pixel 
scale as the EGS (i.e 1.2"/pixel at 24 pm and 4.0"/pixel at 
10pm). In both fields observations reach point source sensitivity 
of ~ 30 pJy and ~ 2 500 pJy at 24 and 70 pm respectively 
(i.e an effective integration time per sky pixel of ~ 35 000s 
and ~ 10 800s respectively). Source extraction was carried out 
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Fig. 1. Photometric accuracy, [S(out)-S(in)]/S(out), derived from the simulation as a function of measured flux density, S(out), for 
the GOODS-N 24 pm (left) and 70 pm observations (right). Shaded gray region shows the distribution that we obtain and the 
solid black lines contain 68 % of the galaxies. Vertical long dashes lines show the limit of our catalog defined as the best balance 
between photometric accuracy, completeness and contamination. Inset plots show the fraction of artificial sources detected in the 
image (completeness) as a function of input flux (orange plain histogram) and the fraction of spurious sources (contamination) as a 
function of flux density (striped black histogram). For the definition of the completeness and contamination see the text. 



using our method with prior positions, leading to a total of 
2151 sources at 24 pm(> 30/iJy, 5<x) and 119 sources at 70 pm 
(> 2.5 mJy) in GOODS-N and 870 sources at 24 pm (> 30yuJy) 
and 44 sources at 70 pm (> 2.5 mJy) in GOODS-S. 



2.1.3. ECDF-S 

The ECDF-S (3 /, 32" ! , -27°48') was observed with the MIPS in- 
strument as part of the FIDEL legacy program (PI: Dickinson) 
and covers a total area of 900 arcmin 2 . In the following we re- 
fer to the ECDF-S field as the outer part of this 30' x 30' area 
while the central part (1 0' x 10') is refered to as the GOODS-S 
field (see Sec tion l2. 1.21 ). Reduced to its outer part, the ECDF-S 
field covers a total area of 520 arcmin 2 (hereafter ECDFS-O). 
The choice to treat separately the GOODS-S and the ECDFS- 
O data is driven by the difference in the 70 pm coverage (the 
central part being deeper than the outer part). The final 70 pm 
mosaic, obtained using MOPEX, has the same pixel scale as the 
EGS and reaches a point source sensitivity of ~ 3.5 mJy (i.e 
an effective integration time per sky pixel of ~ 5 750s). At 24 
pm, the effective integration time per sky pixel being inhomo- 
geneous, the final mosaic reaches point source sensitivity values 
which are function of the source positions. It spans the range of 
~ 40 - 70 pJy over the whole mosaic. 

Using our extraction method we obtained 2185 sources at 
24 pm (> 80//Jy) and 195 sources at 70 pm (> 3.5 mJy). 

2.2. Infrared source detection 

Since most sources in these fields were unresolved, we per- 
formed a PSF fitting technique to extract their photometry. 
Nevertheless, to be able to deal with blending issues and to go 
down to the limit of the data, we used, as a prior information, 
the expected position of the sources. We assume that all sources 
present in the mid-infrared images have already been detected 
at other wavelengths where the depth and the resolution of the 
data are higher. 



For the 24 ^m data, we choose to use the position of the 
IRAC-3.6 pm sources. This choice was motivated by the depth 
of the IRAC-3.6 pm data which is on average a factor of 30 
deeper than MIPS 24 pm data in every field used in this study. 
Hence since we k now that the typic al Supm/S irac ra ti° spans 
the range [3-20] dCharv et all l2004t) . the IRAC data are deep 
enough to contain every 24 pm source. Moreover the resolution 
of the IRAC-3.6 pm data (FWHM ~ 2") being 3 times sharper 
than the MIPS 24 pm data, blended sources separable in the 24 
pm MIPS images (i.e sources separated by more than ~ 2" in 
the 24 pm data) are unblended in the IRAC catalog. 

For the GOODS-N field, we used the publicly available 
IRAC catalog released as part of the GOODS legacy program 
(Dickinson et al., in preparation). This catalog contains 19437 
objects detected at 3.6 pm with a 50% completeness limit of 0.5 
pJy and an astrometric accuracy of 0.37". For the EGS, we used 
the publicly av ailable IRAC EGS c atalog released by the AEGIS 
collaboration (Barmbv et al. 2008). This catalog contains 57,434 
objects detected at 3.6 pm with a 50% completeness limit of 
1 .5 piy. The astrometric accuracy of this catalog is estimated to 
be close to 0.37". Finally, for the GOODS-S/ECDF-S region, 
we used the SIMP LE catalog (Spitzer IR AC / MUSYC Public 
Legacy in ECDF-S; iGawiser et all 2006). This catalog contains 
61233 objects detected at 3.6 pm with a 50% completeness 
limit of 1.5 piy. The astrometric accuracy is the same as in the 
EGS catalog. 

An empirical 24 pm PSF was constructed with isolated 
point like objects present in the mosaic and was then fitted to 
all IRAC positions in the map. The photometry of each 24 pm 
source is defined as the scaled fitted PSF. Then, as for a standard 
aperture correction, the photometry is corrected to account 
for the finite size of the empirical PSF. A residual image was 
created by subtracting from the original image the detected 
24 pm sources. Finally we check in this residual image the 
presence of 24 pm sources missed by the lack of IRAC position 
using DAOPHOT detection. We find no residual sources. 

For the 70 pm data, we choose to start from our previous 24 
pm detections. This choice is straightforward si nce the typical 
Syofim/Si^ttm ratio spans the range [2-100] (Pa povich et alj 
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Field 


Area 


Exposure 


Flux Limit" 


Completeness 11 


Contamination b 


# sources 




(arcmin- ) 


(sec/pixel) 


(ply) 


(% at S = Sib*) 


(% at S = Sitait) 


(with S > S llrait ) 


EGS 24 


548 


14000 


50 (5tr) 


75 


11 


4096 


GOODS-S 24 


91 


35 000 


30 (5o-) 


81 


13 


870 


GOODS-N 24 


194 


35 000 


30 (5cr) 


81 


13 


2151 


ECDFS-O 24 


517 


8 000 


70 (5o-) 


100 





2474 


EGS 70 


548 


7 200 


3000 (60-) 


77 


10 


322 


GOODS-S 70 


91 


10 800 


2500 (6tr) 


73 


12 


44 


GOODS-N 70 


194 


10 800 


2500 (6tr) 


73 


12 


119 


ECDFS-O 70 


517 


5 750 


3500 (6tr) 


85 


7 


195 



a the level of the noise was calculated on the residual images and hence it combines the confusion, instrumental and photon noises. 
b see text for the definition of the completeness and contamination 



120071) and that the 24 fim observations are about 100 times 
deeper. No reliable empirical PSF can be constructed at 70 fim 
since only a few isolated sources can be found in each map. We 
decided to use the appropriate 70 fim Point Response Function 
(PRF ) estimated on the extragalactic First Look Survey (xFLS; 
iFraver et alj l2006ah mosaic and available on the Spitzer web 
site. Then as for the 24 /mi detection, we fit the 70 fim maps 
with this synthetic PSF at each 24 fim position. Once again we 
then check the residual image but no residual sources where 
found, i.e we did not see any evidence for the existence of 70 
fim sources with no 24 fim counterpart at these depths. 

The advantage of this method of priors is principally to 
resolve a large part of the blending issue. Moreover this method 
reduces the cross-matching issue since for each IRAC position 
we directly have a 24 and 70 /mi flux. Then to cross-match our 
mid-infared detections with any other catalog, we can use the 
IRAC astrometry which is more accurate than any 24 or 70 fim 
centroid position. 

To estimate completeness and photometric reliability 
we performed extensive simulations for each field and at each 
wavelength. We added artificial sources in the image with a 
flux distrib ution matching approximately the measure d number 
counts (see IFraver et alj[2006bl : IPapovich et al.l 12004). In order 
to preserve the original statistics of the image (especially the 
crowding properties) the numbers of artificial objects added in 
the image was kept small (as an example for GOODS-N we 
only added 40 sources into the 24 fim image and 4 sources into 
70 fim images). Moreover as the photon noise of the MIPS 24 
and 70 fim data is dominated by the background emission we 
do not need to introduce any additional Poisson noise due to 
the photon statistics of the object itself. We then performed the 
source extraction again to compare the resulting photometry to 
the input values. To increase the statistics we used repeatedly 
the same procedure with different positions in the same field. 
For each field we introduced a total of 20 000 artificial objects. 

In Figure [TJ we show the results of one of the simulations 
obtained for the GOODS-N 24 fim and 70 fim data. We plot 
for each wavelength the photometric accuracy derived from 
the simulation as a function of measured flux density. We 
defined the completeness as the fraction of sources with a flux 
accuracy better than 50%. To estimate the contamination of our 
catalog by spurious sources (i.e sources arising from the noise 
properties of the image and not from real sources) we simulated 
very faint sources which were normally undetectable since they 
were introduced with a flux density fainter than the background 



noise. The contamination was then defined as the fraction of 
these undetectable sources that were extracted by our method. 
Using these simulations we defined the flux density limits of the 
final catalogs as those above which we simultaneously get : 

- a photometric error better than 33% for at least 68% of the 
sources. 

- a completeness of at least ~ 80% (to minimize future 
completeness correction). 

- a contamination by artificial sources lower than 15 %. 
Under this definition, we found that the GOODS-N and 

S 24 fim catalogs, which have about the same exposure time, 
reach a limit of 30 fi]y (see Table 12.2b . We emphasize that 
sources as faint as ~ 15 - 20//Jy can be extracted from the 
24 fim images. However, in order to avoid contamination and 
large incompleteness, both of which are significant systematics 
in the estimation of the LF, we have cut the GOODS 24 fim 
catalogs at a brighter flux density limit. For the EGS 24 fim 
data the flux density limit of the catalog is 50 yuJy. The 24 fim 
exposure map of the outer part of the ECDF-S (i.e ECDFS-O) 
is inhomogeneous, hence we used the lowest exposure time as 
the reference to define the flux density limit. In order to avoid 
inhomogenous completeness corrections we imposed a level for 
the flux density limit such that the completeness is ~ 100%. 
Combining these two conditions, the ECDFS-O 24 fim catalog 
reaches a limit of 70 /iJy. 

For the 70 fim data of GOODS-N and S, the catalog flux 
density limit is 2.5 mJy. For the EGS 70 fim data, the flux limit 
of our catalog is 3 mJy and for the ECDFS-O 70 fim data the 
detection limit of the catalog is 3.5 mJy. 

We would like to emphasize that the calibration factor 
taken to generate the final 24 and 70 fim mosaics are derived 
from stars, whose SED at these wavelengths are generally very 
different from those of distant galaxies. Hence color corrections 
may be necessary to these fluxes (at most ~ 10%). Nevertheless, 
since these color corrections are highly dependent on the 
redshift and the intrinsic SED of the sources, we decided to 
introduce these color corrections directly in the ^-correction 
used to estimate the LF (see Section 13.11 ) since both quantities 
are taken into account in this computation. 



2.3. Spectroscopic redshifts 
2.3.1. EGS 

The optical spectroscopic redshifts of t he EGS field were taken 
from the DEEP2 data release 3 catalog (Dav is et al ] l2007l) . This 
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Field 


Area 


Nb sources 


X-ray AGN 


No X-ray AGN 


# spec-z a 


# phot-z a - b 


# spec-z and/or phot-z a 




(arcmin 2 ) 


24^m/70//m 


24pm/70//m 


24,um/70pm 


24^m/70^m 


24^m/70^m 


24//m/70,um 


EGS 


S48 


4096 H9? 


1 81 111 


^91 S/98S 


1 108/198 

1 1UO / iz,o 




9881 /9S4 








4%/ll% 


96%/89% 


28%/45% 


45%/45% 


73%/90% 


GOODS-S 


91 


870/44 


64/7 


806/37 


371/32 


378/5 


749/37 








7%/16% 


93%/84% 


46%/86% 


46%/ 14% 


92%/ 100% 


GOODS-N 


194 


2151/119 


134/12 


2017/107 


747/72 


1148/29 


1895/101 








6%/10% 


94%/90% 


37%/67% 


57%/27% 


94%/94% 


ECDFS-O 


517 


2474/195 


0/0 


2474/195 


180/16 


1288/168 


1468/184 








0%/0% 


100%/100% 


7%/8% 


52%/86% 


59%/94% 



11 the percentage note in these columns refer to the number of no X-ray AGN 

b number of sources which have a photometric redshift but no spectroscopic redshift 



survey targeted ~50,000 distant galaxies with R < 24.5, using 
the DEIMOS spectrograph on the Keck II telescope. In the 
EGS, unlike the other DEEP2 fields, objects at z < 0.7 were not 
excluded from the selection. The data release 3 catalog contains 
about 31 600 reliable redshifts with 18.5 < R < 24.5. For more 
details on t he sel ection criteria of the DEEP 2 program see 
iDavis et alj d2007l) . 



Q"Az/d +zi = 0.029 with 3.8% of catastrophic err ors dllbert et alj 
2006). As recommended in Ilb ert et all Q2006), we excluded 
all sources with a double peak in their redshift probability 
distribution function (PDF) since most of them did not have a 
reliable redshift. The overlap region between the EGS and the 
CFHTLS D3 reduces the EGS area from 900 arcmin 2 down to 
548 arcmin 2 . 



2.3.2. GOODS-North/South & ECDF-S 

A large number of spectroscopic redshifts has been measured 
for galaxies in the GOODS regions. In GOODS-N a total of 
2376 spectro scopic redshifts come from a combination of vari- 
ous s tudies dCohen etal.ll2000t IWirfh et al] 12004 ICowie et ail 
2004, Stern et al. in prep). In the GOODS-S/ECDF-S region, 
2547 s pectroscopic redshifts were also combined from various 
studie s dVanzella et alj2006tlLe Fevre et al.ll2004lMignoli et alj 
120051) . Although these spectroscopic redshift catalogs were 
obtained from several campaigns with diverse selection criteria, 
when combined together, they provide a rather homogeneous 
sampling. In GOODS-N and GOODS-S, respectively 60% 
and 50% of the galaxies brighter than zab = 23.5 and with a 
photometric redshift lower tha n z p i m t < 1 -2, have a spectroscopic 
redshift (see lElbaz et"aT]|2007l) . 



2.4. Photometric redshifts 
2.4.1. EGS 

In addition to the spectroscopic redshifts we also need 
reliable photometric redshifts. We used publicly available 
photometric redshi fts provided by the TERAPIX and the 
VVDS consortia dllbert et al.l 12006ft using the Canada-Hawai- 
France Telescope Legacy Survey (CFHT-LS). The CFHTLS 
D3 survey consists of a 1 square degree region centered at 
a = 14 ;, 19 m 27*,<5 = +52°40'65" observed through the u, 
g, r, i and z CFHTLS filters. Using this set of photometric 
magnitudes and appropriate spectral energy distributions (SED), 
the "Le Phare" photometric c ode provided a redshift probability 
distribution for each source dllbert et al]l2006l) . 

The final i'ab selected catalog contains 366 030 sources 
with a magnitude limit of i'ab ~ 26.0. For galaxies with 
i'ab < 24 and z < 1.5, these redshifts have an accuracy of 



2.4.2. GOODS-North & South 



For the GOODS fields, we complemented the spectroscopic 
redshifts with photometric redshifts e stimated using Z-PEG 
dLe Borgne & Rocca-Volmerangel 120021) . For both fields we 
used some of the deepest data currently available in the GOODS 
region (UBVRIzJHK, 3.6 /mi, 4.5 /mi) to derive photometric 
redshifts with an accuracy of cr A -/ (1+ ,) = 0.05 - 0.1. For 
more detail on the complete multi-lambda catalog used in the 
computation of the photometric redshifts see Le Borgne et al. 
(In Prep). 



2.4.3. ECDF-S 



For the ECDF-S we used the catalog of photometric red- 
shifts from the COMBO- 17 survey (Class ifying Objects b y 
Medium-Band Observations in 17 filters ; IWolf et alj 12004). 
The COMBO- 17 survey observed a 31.5 x 30 arcmin 2 area 
with 5 broadband and 12 narrowband filters. The catalog 
contains 63 501 objects selected on the deep R-band image 
(R < 26) carried out with the Wide Field Imager (WFI) at the 
MPG/ESO 2.2-m telescope in La Silla. Using an extensive SED 
library, they fitted the photometry of each object to compute its 
redshift probability distribution. Finally the COMBO-17 catalog 
contains 24 216 objects with a reliable photometric redshift. 

For sources with R < 24 and z < 1.2, these redshifts have 
an accuracy of cr^ z /^ +Z ) = 0.1. For sources fainter than R = 24 
the redshift accuracy drops dramatically, therefore we did not 
consider in the following such faint objects. 
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2.5. Removing Active galaxies 

The infrared luminosity of a normal star-forming galaxy can be 
severely affected by the presence of an Active Galactic Nucleus 
(AGN). Hence to study the star formation history of normal 
galaxies it is necessary to identify and remove these AGNs. 

In order to test the presence of AGNs in the GOODS 
field, we used the deepest Chandra X-ray observations i.e the 
lMs maps for the CDF -S and the 2 Ms maps for the CDF-N 
dAlexander et al.l 120031) . For both fields we used the corre- 
sponding public X-ray catalog. For the EGS field we also used 
Chandra observations taken from the public AEGIS-X data 
release 2 point so urce catalog which has a 200 Ksec coverage 
dLairdet al.ll2008h . Finally for the ECDF-S we used the 250 
Ksec Chand ra observations, wh ich flank the 1 Ms CDF-S 
observations (Lehmer et al. 2 0051) . 

AGNs were identified as galaxies with either 
L z [0.5 - 8.0 keV] > 3 x 10 42 ergs" 1 and/or a hardness 
ratio greater than 0.8 (ratio of the counts in the 2-8 keV 
to 0.5 - 2 keV passb ands) dBauer et"aT. 1 12004b . As mentioned 
in Ba uer et a D (12004 . some AGNs might be missed by this 
technique. The fraction of such galaxies remains a matter of 
debate until today but the global consistency of the 24-70 fim 
color of the galaxies studied here (as discussed in Section I4.21 > 
suggests that they should not play a dominant role (see also 
Fadd a et al.1l2002l) . We did not subtract the AGN contribution 
to the IR light of those galaxies harboring an AGN, instead we 
conservatively assumed that any galaxy contaminated by an 
AGN was not forming stars. It is well-known that AGNs do also 
harbor star formation, especially those emitting strongly at IR 
wavelengths, but such subtraction would be highly speculative 
at the present level of our knowledge. 



2.6. The final IR galaxy sample 

To construct our final catalog, we first truncated our mid- 
infrared catalog to the region covered by both Spitzer and 
the optical redshift surveys used in this study. This step was 
especially important for the EGS field where the total FIDEL 
area of 900 arcmin 2 is reduced to 548 arcmin 2 to overlap the 
CFHT-LS survey (i.e the photometric redshift survey). We then 
cross-matched the 24 fim catalog with the X-ray observations 
to exclude AGNs. We found that ~ 6% (~ 10%) of the 24 fim 
(70 fim) sources were affected by an X-ray AGN and hence 
removed from our catalog. The rest of the data was then cross- 
matched with the different spectroscopic catalogs and finally the 
remaining sources were matched with the photometric catalogs. 
All these cross-matching steps were made using the IRAC 
astrometry of the mid-infrared sample and a matching radius 
of 1.5 arcsec. The choice of 1.5 arcsec is motivated by the size 
of the IRAC-3.6jum PSF FWHM (~ 2"). In case of multiple 
associations we selected the closest optical counterparts to the 
IRAC centroid. We found that 8%, 10% and 9% of the MIPS 
24 fim sources have multiple associations with optical sources 
brighter than I < 26 mag for the EGS, GOODS-N, GOODS-S. 
For the ECDFS-O we found that 1% of sources have a multiple 
association with optical sources brigther than R < 24 mag. 

All the different steps described previously and applied to 
the different fields of our study are listed in the Table [2] The 
redshift completeness is larger than 90 % (70 % except for the 
ECDFS-O) at 10 fim (24 fim) for the sources with no AGN 
contribution. 
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Fig. 2. Comparison between photometric and spectroscopic red- 
shift for 2406 galaxies of our 24 yum-selected catalog. Dashed 
lines represent the relative errors (i.e a-^zi(\^z) = 0.06). 77 gives 
the fraction of source with Az/(1 + z) > 0.15. 



2. 7. Redshift uncertainties and sample completeness 

One of the major concern of all studies heavily relying on 
photometric redshift is their possible redshift incompleteness 
for distant and/or faint galaxies. As these issues may have a 
nonnegligible effect when estimating luminosity functions as a 
function of redshift, one should pay attention to the consequence 
of these limitations. 

In Section 12.41 we have quoted the accuracy of all photo- 
metric redshift catalogs used in this work. Nevertheless since 
these accuracies were based on an optically-selected sample 
we need to characterize the typical uncertainties of our 24 
yizm-selected sample. In Figure [2] we compare the spectroscopic 
and photometric redshifts of all 24 fim sources with a spectro- 
scopic redshift after combining the four fields (2406 sources, 
see Table We observe a very good agreement between the 
photometric and spectroscopic redshifts with an accuracy of 
o"Az/(i+z) = 0.06 and a median value of -0.005. We notice that 
this accuracy is not statistically larger than the accuracy quoted 
for the optically-selected sample. We also found 4% of catas- 
trophic objects defined as those objects with Az/(1 + z) > 0.15. 
Even if those galaxies are statistically rare, they could have a 
major effect on the estimates of the LF by boosting the bright 
end slope. Therefore in the following study we will pay close 
attention to quantify their effect on the estimates of the total 
infrared LF (see Section [3~TT i. 

As previously mentioned and quoted in Table [2] a high 
fraction (~ 80%) of our final IR galaxy sample have been 
associated with a spectroscopic and/or a photometric redshift. 
Nevertheless the remaining fraction of sources with no redshift 
association are not randomly distributed in the redshift space. 
As a result, we need to quantify the redshift above which 
the photometric redshift incompleteness affects our study. In 
Figure [3] we present the photometric redshift completeness as a 
function of the optical magnitude (solid lines). These histograms 
allow us to estimate the optical magnitude above which the 
determination of a reliable photometric redshift start to be 
problematic. Independently, we also present in figures [3] the 
optical magnitude distribution of the 24 fim sources for different 
redshift bins. For clarity these optical magnitude distributions 
have been scaled with an arbitrary factor. By combining these 
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Fig. 3. Fraction of 24 /an sources with an optical counterpart 
which have a reliable photometric redshift as function of optical 
magnitude (solid line). The dashed line, in the EGS panel, rep- 
resents the redshift determination completeness as a function of 
the optical magnitude before the exclusion of all sources with a 
double peak in their PDF (see Section l2~4l i. The shaded regions 
show the optical magnitude distribution for various redshift bins 
(all histograms have been scaled with an arbitrary factor). Note 
that the ECDFS-0 panel uses the R-band magnitude. 



two informations (i.e the redshift determination completeness 
as a function of the optical magnitude and optical magnitude 
distribution) one can have an estimates of the redshift incom- 
pleteness of his sample. For example, if the photometric redshift 
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Fig. 4. Redshift distributions of 24 jjm sources selected with 
F 2 4fim > 70 //Jy in the GOODS-S and N fields (solid line), in the 
EGS field (dash line) and in the ECDFS-O field (dot line). 



catalog is highly incomplete for sources fainter than R > 24 and 
in the same time a large fraction of the 24 fim sources that lie 
at 1 < z < 1.3 have R ~ 24.5, one can expect a large redshift 

incompleteness in this redshi ft bin. 

As already discussed in lLe Floc'h et a I] (120051) we observe 
that sources situated at higher redshift have fainter optical 
counterpart and that the redshift determination completeness 
decreases at faint magnitudes. For the GOODS fields we notice 
that the depth of optical observations and the use of deep near 
infrared data allow us to have a complete redshift association up 
to z ~ 2. Indeed we notice that the redshift completeness of the 
GOODS field falls at optical magnitudes (i.e R > 27.5) where 
the distribution of 24 fim sources that lie at 1.3 < z < 2.5 is 
already low. Therefore for the GOODS fields, sources with no 
redshift association (~ 5%) are likely located at z > 2 and hence 
should not affect the study of the evolution of the total infrared 
LF up to z = 1.3. For the EGS the redshift determination 
completeness highly decreases at an optical magnitude range 
(i.e I-mag~ 22.5) where a large fraction of sources are located 
at z < 1.3. This trend can lead to a high incompleteness of 
sources located at 1 < z < 1.3. However we notice that the 
decrease of this redshift determination is driven by sources with 
a double peak in their probability distribution functions (PDF; 
see the dash lines in Figure |3). All galaxies for which we have 
both a double peak distribution in their photometric redshift and 
a spectroscopic redshift are located at z > 1 . The minor effect 
of "double peak" galaxies on the redshift distribution can be 
inferred by the comparison of the redshift distribution of the 
EGS and the GOODS fields (for which we expect no redshift 
incompleteness up to z ~ 2; see Figure [4] and Figure |3). Indeed 
we do not observe large differences in the redshift distribution 
of these two fields between < z < 1.3, whereas at z > 1.3 
the number of sources in the EGS drops to zero likely due to 
redshift incompletene ss. Finally for the ECD FS-O field we find 
the same results than iLe Floc'h et alj d2005l) . i.e the magnitude 
distribution of the two highest redshift bins (i.e 0.7 < z < 1.0 
and 1.0 < z < 1-3) are truncated at the faint end which 
indicate that these redshift bins are affected by incompleteness. 
Nevertheless we can estimate from the shape of the magnitude 
distribution that only ~ 15% of sources should be affected by 
this incompleteness and hence no large effects are expected 
on the estimates of the infrared LFs. Indeed we will see in the 
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following that the infrared LFs inferred from the ECDFS-O and 
GOODS fields are, within the error bars, in very good agreement 
(see Section [3] |4] and |5). To conclude we infer that the redshift 
incompleteness should not affect the conclusions of our study 
since the uncertainties corresponding to this misidentification 
are lower than the uncertainties introduced by the ^-correction 
and the photometric redshift measurements. 

Assuming that our photometric redshift catalogs do not 
suffer from incompleteness up to z = 1.3, we can construct a 
reference sample of galaxies located at z < 1.3. When compared 
to this reference sample, we find that the fraction of sources 
with a spectroscopic redshift is 56%, 35% and, 9% for the 
GOODS, the EGS and the ECDFS-O fields respectively. 



3. The rest-frame 15 Micron Luminosity Functions 

In this Section, we present the redshift evolution of the rest- 
frame 15 /mi luminosity function (LF) that we derive from the 
catalogs described in the previous section using the observed 
MIPS 24 /mi flux density. The 15 /mi LF is computed for three 
redshift bins centered at z=0.55, 0.85 and 1.15 (same bins as in 
Figure for which the 24 /mi observations correspond to the 
rest-frame 15, 13 and 11 /mi wavelengths respectively. The ex- 
trapolation that needs to be applied to compute the 15 /mi lumi- 
nosity is therefore negligible in the first two redshift bins (still 
it is carefully computed here) and becomes larger in the high- 
est redshift bin due to the fact that the observed 24 /mi is more 
distant from the rest-frame 15 yum. While the derivation of the 
rest-frame 15 /mi LF has already be en addressed in some pre- 
vious studies dLe Floc'h et al. 1 12005b . this is done here for two 
reasons. 

First, we wish to check and demonstrate the consistency of 
our technique and catalogs with those previous studies at a com- 
mon wavelength before deriving the rest-frame 35 /mi LF. As 
we show it here, our results are not only consistent with those 
previous studies but even extend them to fainter luminosities as 
well as they provide a more robust constraint on the 15 /mi LF 
because ou r catalogs are three tim es deeper at 24 /mi than the 
one used in Le Flo c'h et all d2005l) . We are also able to validate 
the results against the effects of cosmic variance by comparing 
three independent fields: the EGS, the GOODS and the ECDFS- 
O fields. 

Second, this study will serve as a basis for the extension of 
the derivation of the faint end of the 35 /mi LF. For that purpose, 
we will first determine the typical ratio of the observed 70 over 
24 /mi flux densities per redshift and flux density bins in order 
to associate 24 //m sources to a given flux density bin at 70 /mi. 
Then we will derive the contribution of this bin to the 35 /mi LF. 

The sizes of the redshift bins [0.4-0.7], [0.7-1.0], [1.0-1.3] 
were defined such as to encompass a large enough number of 
sources in the computation of the 35 /mi LF. For consistency, we 
used the same redshift bins here for the 15 //m LF. 



3.1. Methodology 

We define here the k-correction as the ratio between the rest- 
frame 15 and 24/(1 + z) /mi luminosities as in Eq.Q] 

vi5/j m ii 5/ "" = v 2 4 tim^-ndlS (24fim)k(z) (1) 

where L l v 5i " n is the rest-frame monochromatic luminosity in 
WHz S(24 jt/m) is the observed 24 /mi flux density in 




0.6 I , , , i , , , i , , , i , , , i , I 

0.4 0.6 0.8 1.0 1.2 

Z 

Fig. 5. 15 /mi ^-correction as a function of the redshift (as in de- 
fined in Eq.Q~|l for the MIPS-24 /mi pass band. The thin and the 
thi ck solid lines corresp ond to the ^-correction obtained using 
the lCharv & ElbazM2001l) 10 11 L Q and 10 12 L G templates respec- 
tively. The thin and the thick dashed lines corres pond to the k- 
correction obtained using the lLagache et al.l (120031) 10 11 L Q and 
10 12 L Q template respectively. Finally the thin and thick dashed- 
dot-dot lines corresp ond to the ^-correction obtained using the 
iDale & Helov3 (120021) 10 11 L Q and 10 12 L Q template respectively. 

Wm~ 2 Hz _1 , di is the luminosity distance and k(z) is the k 
correction. This correction depends on the exact spectral energy 
distribution (SED) of the galaxies. In the following, we will 
use the Chary & Elbaz (2001; hereafter CE01) library of 105 
template SEDs as proxies of the typical SED of galaxies as- 
suming that these SEDs, constrained by local galaxy properties, 
remain valid in the redshift range considered here. There are 
good reasons to believe that this is indeed the case, such as the 
consistency between the total IR luminosity derived from the 
observed mid-IR and radio flux densities of galaxies (Elbaz 
et al. 2002, Appleton et al. 2005, Marcillac et al. 2006). Until 
enough data are gathered combining broadband photometry and 
spectroscopy from IRS (the Spitzer infrared spectrograph), we 
use three available libraries of template SEDs to quantify the 
level of variations of the k-correction in this redshift range at 
first order. These three librarie s are the CE01, LDP (Lagache, 
Dole & Puget 2003) and DH dDale & Heloul l2002t) ones. We 
convolve the redshifted SED templates with the Spitzer 24 /mi 
filter to compute the corresponding ^-corrections. In Figure 
[5] we present the inferred ^-corrections for the three different 
libraries. We notice that for a given total infrared luminosity 
(i.e L/r = L(8 - 1000 /mi)), all three libraries give a slightly 
different ^-correction as a function of redshift. Nevertheless 
these discrepancies reach at most a value of 20%. In Table [3] we 
present the average value of the ^-correction inferred for our 
sample in the three redshift bins for the three different libraries. 
We notice that their average values are of the same order with 
a typical offset of 10%, 5% and 5% in the three redshift bins 
respectively. Hence, the choice of a particular SED library 
should not affect much the computation of the rest- frame 15 
microns LF. 

In the following, the 15 microns luminosity (vL v /Jm ) of each 
galaxy is derived using the CE01 library. This choice was made 
to be homogeneous with the rest of the study since we will also 
use the CE01 templates to derive the total infrared luminosity of 
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0.4 < z < 0.7 


0.7 < z < 1.0 


1.0<z< 1.3 




CE01 


LDP 


DH 


CE01 


LDP 


DH 


CE01 


LDP DH 


Median 


1.29 


1.16 


1.16 


0.88 


0.90 


0.89 


0.98 


0.95 0.97 


Mean 


1.32 


1.15 


1.15 


0.87 


0.91 


0.90 


1.00 


0.97 0.99 


Sigma 


0.20 


0.06 


0.07 


0.11 


0.08 


0.07 


0.15 


0.10 0.11 



galaxies. Indeed, as illustrated in Section l4~2l the CE01 library 
appears to best reproduce the observed luminosity-luminosity 

correlations even at high redshift (i.e z ~ 1). 

T he LF are derived using the 1/V max method (Schmidt 
1968) which does not require any assumption on the shape 
of the LF. The LF is then directly derived from the data. The 
comoving volume associated to any source of a given luminosity 
is defined as V max = V zmax - V zm ,„ where zmin is the lower 
limit of the redshift bin, and zmax is the maximum redshift at 
which the object could be seen given the flux density limit of 
the sample, with a maximum value corresponding to the upper 
redshift of the redshift bin. For each luminosity bin, the LF is 
then given by 

^ L ( . v max,i 

where V max is the comoving volume over which the z'th galaxy 
could be observed, AL is the size of the luminosity bin, and w, 
is the completeness correction factor of the ith galaxy, w, equals 
1 for sources brighter than S 24^,1 ~ 100/iJy and decreases at 
fainter flux densities due to the incompleteness of the 24 fim 
catalog (see the red histogram of the inset plot in Figure [TJ. 
These completeness correction factors are robustly determined 
since they only affect the faintest luminosity bin and by at most 
a factor of 0.7. Hence none of the conclusions presented here 
strongly rely on this correction. 

The four different fields studied here (EGS, GOODS-N, 
GOODS-S and ECDFS-O) have different flux density lim- 
its, therefore we will treat separately the GOODS-N and S 
(S'™« m = 30 /Jy), the EGS (S'>™' m = 50juJy) and the ECDFS-O 

0^24#m = 70 /Jy). This will allow us to test the role of cosmic 
variance. 

The uncertainties of the LF values depend on the numbers 
of sources in each luminosity bin, on the photometric redshifts 
uncertainties, the ^-correction uncertainties and finally on the 
24 and 70 //m flux uncertainties. To compute these errors we 
performed Monte Carlo simulations by creating 1000 simulated 
catalogs. These simulated catalogs contain the same number 
of sources as the original one but we attributed to each source 
with a photometric redshift a new redshift randomly selected 
within its probability distribution function. Sources with a 
spectroscopic redshift have been left unchanged. To take into 
account the redshift catastrophic errors we randomly selected 
4% of sources with a photometric redshift (corresponding to 
the fraction of catastrophic sources found in Section [2T7T > and 
then attributed to these sources a random redshift selected 
assuming a flat distribution between < z < 3. We notice that 
assuming a flat distribution may not be totally realistic since 
catastrophic errors are generally bimodal (e.g, assigning to a 
low redshift galaxy a much higher redshift due to mistaking 
one spectral break for another). Nevertheless the sample that 
we use to compare the photometric and spectroscopic redshifts 
is not large enough to describe rigorously this bimodality. 



However, assuming a flat distribution may be a more conser- 
vative assumption since this can yield an overestimation of the 
uncertainties due to redshift catastrophic errors. Indeed if the 
bimodal distribution of the catastrophic errors is characterized 
by Az > 3 (i.e due to mistaking the Lyman break for a Balmer 
break) then it will only affect the estimates of the LF at z > 3. 
To take into account flux uncertainties we attributed to each 
source a new flux selected into a Gaussian distribution centered 
at the original flux of the source and with a dispersion computed 

as CT = ^((T calibration) 2 + {& source) 1 where CT C atib r ation is the 

calibration uncertainty estimated to be ~ 10% of the source 
flux density and o- source is the source photometric uncertainty 
computed on the residual image. Finally to take into account 
the ^-correction uncertainties, we selected the source luminosity 
assuming a Gaussian distribution centered at the original source 
luminosity and with a dispersion of 0.05 dex. Using these 1000 
simulated catalogs, we then computed again the rest frame 15 
yizm LF and defined the total uncertainty in each luminosity bin 
as the quadratic sum of the Poissonian error (oc VI INb sources ) 
and the dispersion given by the Monte-Carlo simulations. We 
found that the Poisson statistics dominate the total uncertainties 
at the bright luminosity bin whereas at faint luminosities the 
total uncertainty is dominated by the Monte Carlo uncertainty. 

To study the evolution of the rest-frame 15 microns LF 
between z ~ and z ~ 1.3 we u sed as a local reference the 
local 15 microns LF of [Xu! (2000) derived using the bivariate 
(15 yum vs. 60 /im luminosity) method and based on ISOCAM 
observations. This LF is then fitted with a double power law 
function. This choice is driven by the fact that the ultimate 
goal of this paper is to derive the evolution of the total in- 
frared LF which is itself bes t fitted with a double power law 
function (ISanders et al.l 120031) . Hence to be homogeneous we 
have made the choice, in the following, to fit the 15, 35 fim 
and total inf r ared L Fs with double power laws. Using the 
ISanders et all d2003l) best-fit (i.e </> oc L ° 6 for log(L/L )< 10.5 
and (p K Zr 215 for log(L/L )> 10.5) and a bivariate method, 
we derived the values of the two power laws for the 15 /im LF. 
We find that <p oc L~ 5 for log(L/L )< L k „ ee and cp oc L 2 15 for 
log ( L/L ) > Lfa, ce . Then, using a% 2 minimization, we fitted the 
Xu (2000) observations with this function, leaving Ly, ee and 
the normalization as free parameters. We also fit independently, 
with the same method and the same function, the rest-frame 15 
microns LF derived in the three redshift bins (see Tabled. 

In Figure [6] we present the rest-frame 15 fim LF derived 
with the \ IV ma X method for the three redshift bins considered 
in this study. As already mentioned, these LFs have been 
derived jointly for the GOODS-N and S fields and separately 
for the EGS and ECDFS-O fields. The negligible role of cosmic 
variance is illustrated by the very good agreement between the 
three solutions in the derived 15 /im LF. 
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Fig. 6. Rest-frame 15 microns LF estimated for three redshift bins with the l/V„ mx method. The light red squares are obtained for 
the combined GOODS-North and South fields. The dark blue squares an d the light blue s quare s are obtained for the EGS and the 
ECDFS-O fields respectively. Empty triangles show the LFs der ived bvlLe Floc'h et al.l ((2005) in the ECDFS at redshifts ~ 0.7, 
~ 0.9 and ~ 1.1. Asterisks show the local reference taken from IXul (120001) and the dotted line is a fit to this data with a double 
power law function (see text). In the panel showing the 0.7 < z < 1.0 redshift bin, we reproduced the best fit of the LF obtained 
in the two other redshift bins (i.e 0.4 < z < 0.7 dot-dashed line and 1.0 < z < 1.3 triple-dots-dashed line). The shade area in each 
panel spans all the parametric solutions obtained with the x 2 minimization method and compatible with the data within 1 <x. The 
inset plot represents the evolution of the <pk„ ee and Lk m e as function of redshift obtained from fits to all field {red fill circle), GOODS 
(empty star), EGS (empty triangle) and ECDFS-O (empty square). The redshift evolution of L^ nee and <pknee was fitted using laws in 
the forms of (1 + z)° L and (1 + z) a *, respectively. 



3.2. Results 

The resulting rest- frame 15 pm LFs are presented in Figure [6] 
They provide three significant results. 

First, as written in the previous section, no strong effect of 
cosmic variance (i.e ~ 0.05 - 0.1 dex) is found between the 
three fields, which provide the same 15 pm LF within the error 
bars. This observed estimate of the cosmic variance is totally 
consistent with theor etical estimates of ~ 0. 07 dex based on the 
area of our fields (see lTrenti & Stiavellill2008h . 

Second, the comparison of the rest-fra me 15 LF 
obtained in this study and the ones inferred by iLe Floc'h et al.l 



(2005) are in agreement for all 3 redshift bins. We note that 
the inferred LF was here extended to fainter luminosities using 
our ~ 3 times deeper catalog which allowed us to obtain better 
constraints on the faint end slope of the LF. 

Third, we find a substantial evolution of the LF with 
redshift. This evolution, already revealed by previous studies 
(ILe Floc'h et al.ll2005t ICharv & Elbazll200Tt ICaputi et al.|[ 2007 ; 
IXul 120001: lLagache et al.l 120031: iFranceschini et al.l 1200 ll) is 
usually expressed through an evolution both in density and 
in luminosity of the local LF (i.e p(L,Q)). Assuming that the 
shape of the LF remains the same since z ~ 0, these studies 
express the evolution as p(L,z) = g(z)p(L/f(z),0), where g(z) 
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Fig. 7. 35 /mi ^-correction as a function of the redshift for the 
MIPS-70 /mi pass band. Lines are as in Figure[5] 



and f(z) describe the density and the luminosity evolution 
through g(z) = (1 + z) p and f(z) = (1 + z) q . The weakness of 
this parameterization is that it assumes a monotonic evolution 
of the LF preventing any turn-over in the evolution (e.g an 
increase in density that would be followed by a decrease). 
This kind of tur n-over has already been seen at other infrared 
wavelengths (e.g Caputi et al. 2007). Therefore, we chose to fit 
independently the 3 different redshift bins with double power 
laws. Then tracking the value of Ltme and <pknee (see Sect. 13. It 
as a function of redshift enabled us to constrain the evolution 
of the LF. We find that the variation of <pknee with redshift is 
not large and starts with a small rise from z ~ to z ~ 0.55 
followed by a small fall up to z ~ 1.15. Hence the 15 /mi LF 
appears to be rather consistent with the standard definition of 
a pure, monotonic, luminosity evolution in this redshift range, 
proportional to f(z) = (1 + z) iJ±03 . We note that the variations 
in density or of ^hwe are small and of the same order than the 
difference in the galaxy number densities found between the 
GOODS and EGS fields in the faint end of the LF. Such small 
variations could well be produced by cosmic variance. 



4. The rest-frame 35 Microns Luminosity Functions 

We now compute the rest-frame 35 /mi LF, based on 70 /mi 
observations, to minimize ^-correction in the z < 1 .2 range and 
hence reduce the dependency of the LF on the assumed SED 
library. 



4.1. Methodology 

In order to minimize ^-correction uncertainties due to po- 
tentially evolving SED with redshift, we chose to derive the 
FIR LF at the closest rest-frame wavelength associated to the 
MIPS 70 /mi filter in the two most distant redshift bins (i.e 
0.7 < z < 1 .0, 1 .0 < z < 1 .3), which translate into the rest-frame 
30-40 /mi wavelength range. Therefore we chose to derive the 
rest-frame LF at 35 microns. In the lowest redshift bin (i.e 
0.4 < z < 0.7), the MIPS 70 //m filter probes the rest-frame 
40-50 /mi wavelength range where the ^-correction is larger but 
local SEDs such as those u sed here were prove d to be quite 
robust at these low redshifts fflavouzet et al.ll2008h . 

As already discussed in this paper, the three different SED 
libraries give ^-corrections which agree within 10% on average 
and up to 20% at most around z ~ 0.5 (see Figure [7] and Table 
HJ. In the following we use the CE01 library and we will show 
in Section [4~2l that this is indeed the best one, among the three 
considered here, for our purpose. 

The rest- frame 35 /mi LF was derived using the 1/V max 
method and the fc-corrected 70 /mi data. The completeness 
correction factors were derived using the simulations described 
in Section [2T2J The GOODS-S and N fields are treated jointly 
(flux limit ~ 2.5 mJy), the EGS and the ECDFS-O were treated 
separately (flux limit ~ 3 mJy and ~ 3.5 mJy). As for the 15 
microns LF, to take into account most of the uncertainties, we 
performed Monte Carlo simulations. For these Monte Carlo 
simulations the 35 //m ^-correction uncertainties were assumed 
to be of 0.08 dex. This conservative value was derived using the 
largest disagreement observed between the three SED libraries. 
Indeed, as shown in Figure[7] we observe at z 0.6 a disagreement 
of 0.08 dex between the LDP and DH libraries. 

There is no direct local reference available in the literature 
for the 35 /mi LF. Therefore we construct ed a local re f erence 
from the rest-frame 25 microns LF of IShupe et alj (QUI 
derived using the IRAS Faint Source Survey. Using a bivariate 
method and the CE01 library, we converted the 25 microns 
LF into a rest-frame 35 microns LF. As already mentioned, 
in the sake of consistency all along this study, we fitted this 
local reference and the three other LFs (corresponding to 
our three redshift bins) with double power laws with fixed 
slopes, i.e <p oc Zr - 6 for log(L/L )< Lk nee and <p oc Zr 2 - 2 for 
log(L/L )> L knee (see Tabled 

We have computed the rest-frame 35 /mi LF in two suc- 
cessive steps. As discussed above, we started with the same 
technique as for the rest-frame 15 /mi LF, using the catalog of 
sources detected at 70 /mi and correcting for incompleteness, 
k-correction and 1/Vmax. The results are presented in the 
Figure with filled stars for the EGS, GOODS (N & S) and 
ECDFS-O fields. For each redshift bin, we fitted the rest-frame 
35 //m LF using the local double power law reference with fixed 
slopes. In Figure [9] the light shaded areas represent the areas 
spanned by all the parametric solutions which are compatible 
within 1 <x with the rest-frame 35 //m derived using the 70 //m 
observations. Here again, we find that the bright end of the 
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LF is well constrained and does not show any evidence for a 
change of slope. The regime that we were able to probe down 
to the 70 fim detection limit of 2 mJy corresponds to total IR 
luminosities of 0.7, 2 and 5x10' 1 Lq computed from the CE01 
library. Hence, apart from the upper redshift bin, this luminosity 
domain covers the full range of LIRGs and ULIRGs. This brings 
direct constraints, at these redshifts, to the relative contributions 
of LIRGs and ULIRGs to the comoving SFR density. Previous 
studies, based either on extrapolations from the mid-IR or on 
uncertain dust extinction corrections from the UV, were likely 
less robust than the derivation of a total IR luminosity from 
rest-frame 35 fim luminosities. In a second step, presented in the 
next section, we extend our derivation of the 35 fim LF below 
its knee by using a stacking analysis. 



4.2. Constraining the faint end 35 microns LF through 
stacking 

In order to stack positions of 35 fim sources below the detec- 
tion limit, it is necessary to adopt some prior positions. The 
24 fim catalogs provide the best proxy for this purpose since 
the 24 fim images are two orders of magnitudes deeper than the 
35 fim ones and the two wavelengths are the closest available. 
If the IR SED of galaxies of a given luminosity would not vary 
much since z ~1, we would expect a direct correlation between 
the 24yum/(l+z) and 70 fim/(l+z) luminosities since local galax- 
ies do follow tight correlations between their 6.75, 12, 15, 25, 
60 and 100 fim luminosities (see e.g. Chary & Elbaz 2001). We 
have tested the existence of such correlations both in the bright 
luminosity range where galaxies are detected at 24 and 70 fim 
and in the lower luminosity range by stacking 70 fim images on 
the positions of sources detected only at 24 fim. 

We started by dividing the 24 fim catalogs into three sub- 
samples, one for each redshift bin. Then, we separated 24 fim 
sources per luminosity bins of 0.4 dex. For each 24 fim luminos- 
ity bin we stacked on the residual 70 fim image all sources with 
no 70 fim detection. The photometry of the stacked image was 
then measured using an aperture of 16", a background within an- 
nuli of 40 and 60"and an aperture correction factor of 1.705 (as 
discussed in the Spitzer observer's manual). Finally, the mean 
70 fim flux density (F? 0/o ") for a given 24 fim bin was computed 



following Equation [3] 



^70 fim _ 
bin 



m x f™7 + y" , f: 

slack ^i~Y i 



70/JH7 



n + m 



(3) 



where F 1 ^? is the 70 fim flux density, obtained through 
stacking, of all the 24 fim sources undetected at 70 fim (sample 
which contains m sources) and F ■ is the 70 fim flux density 
of the /th 24 fim sources detected at 70 fim (sample which 
contains n sources). This procedure was performed using a 
sliding 24 fim luminosity bin and a step of 0. 1 dex. The result 
of this stacking analysis is shown with filled red diamonds in 
Figure[8] We notice that for the brightest 24 fim luminosity bins, 



bin 



is totally dominated by the flux of detected sources since 
only a small fraction of the sources are not detected at 70 fim. 
On contrary for the faintest 24 fim luminosity bins, F 701 "" is 

dominated by F 1 ^. 

Uncertainties of stacking values depend on the homogeneity 
of the underlying sample and on background fluctuations. To 
estimate the uncertainties related to the sample inhomogeneity 
(i.e (Tinhomo) we perform a standard bootstrap analysis. To 



estimate the uncertainties due to background fluctuations (i.e 
Chuck) we perform a stacking analysis at random positions using 
the same numbers of stacked elements. This stacking analysis 
at random positions is repeated 100 times and then cr iflci . is 
given by the dispersion of these stacking fluxes. Finally the total 
uncertainty is defined as the quadratic sum of <Ti„homo and (Tback- 
These uncertainties are shown in Figures [8] 

The first result that we find here is that in both regimes 
- i.e. where sources are strong enough to be detected at both 
wavelengths and where sources are only seen at 24 fim - the 
two luminosities in the rest-frame wavelengths 24jum/(l+z) 
and 70 fim/(l+z) do follow a correlation. This correlation is in 
agreement with the CE01 library but is a factor of two below 
predictions from the LDP and DH libraries. This correlation 
being found in the three redshift bins it confirms that one can 
complement the 70 fim LF at the faint end based on 24 fim 
detections. Therefore using the inferred 24/70 fim correlation 
(red line in Figure [8]l we have derived the 70 fim luminosity 
of each source starting from their 24 fim flux density and then 
produced the full 35 fim LFs. These LFs are represented with 
filled squares in Figure [9] Using the local double power law 
reference with fixed slopes, we then fitted these rest-frame 
35 fim LFs and all solutions of this fit compatible, within 1 cr, 
with the data are presented as the dark shaded area of Figure 
[9] As revealed by the overlapping of the light and dark shaded 
areas, the resulting 35 fim LF are fully consistent with the ones 
based on direct 70 fim detections and extend them far below the 
knee of the LF. 



4.3. Results 

We used the combination of direct 70 fim detections and 24 fim- 
derived luminosities, converted into 35 fim luminosities from 
stacked 70 fim images, to derive 35 fim LF over more than two 
orders of magnitude, allowing us to position the knee of the LF 
in all three redshift bins (see Figure [5). The data do not show 
any obvious sign for a change of slope of the LF both at the 
faint and bright ends, from z ~0 to 1.3. This confirms that our 
choice to keep the same double power-law slopes than the ones 
used locally is pertinent. Like for the 15 fim LF, we only find 
marginal evidence for a variation of tpknee with redshift, while 
Lknee follows a rapid evolution proportional to (1 + z) 3 • 6±0 - 5 law 
from z ~ to z ~ 1.3, consistent with pure luminosity evolution. 
The flat slope of the bright end of the LF, typical of the IR LF, 
and its rapid evolution with redshift clearly illustrates the rising 
power of luminous galaxies over less luminous ones when going 
to higher redshifts. In the following section, we discuss the im- 
plications of the redshift evolution of the 35 fim LF on the total 
IR LF. 



5. Total Infrared Luminosity Function 

In this section, we use the best proxy to compute the total IR 
luminosity of galaxies (Lir), i.e. the observed 70 fim data when 
available or else the 24 fim data. Lir is computed over the wave- 
length range 8-1000/im. Under this definition, the SFR and 
Lir are related by Eq. [4] (Kennicutt 1998) for a Salpeter IMF 
(0(m)oc irr 2 - 35 , between 0.1-100 M ). Note that we integrate 
the SED curve to derive Lir instead of using the classical ap- 
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Fig. 8. The 70 Vs. 24 fim correlations as revealed by the observations and our stacking analysis in the three redshift bins considered 
in this study. The empty and filled square present the 70 Vs. 24 yum correlations observed for sources which have a photometric 
and a spectroscopic redshift respectively. The red diamond show the results obtained using our stacking analysis (see text). The 
solid lines, the dashed lines and the triple-dots-dash lines represent the expected correlations for the CE01, the LDP and the DH 
library, respectively, at the lowest and the highest redshift of each redshift bin. The red solid line represents the inferred 24/70 //m 
correlation. At the bottom of each plot we present the fraction of 24 yum sources detected at 70 yum as a function of the 24/(1 +z) yum 
luminosity. 



proximation using the fou r IRAS broadband filters as defined in 
ISanders & Mirabell d 19961) . 

SFR [Mq yr 1 ] = 1.72 X 10- 10 L IR [L Q ] (4) 

Up to a redshift of z ~ 1, the local vL v -Lir correlations appear 
to remain valid (see Bavouzet et al. 2008, Elbaz et al. 2002, 
Appleton et al. 2004) and the wavelengths closest to the peak 
far IR emission provide the most accurate estimators of Lir 
(Bavouzet et al. 2008). The k-correction factors required to 
convert a measurement at 70 yum into an Lir present a smooth 
evolution with redshift, consistently followed by all three SED 
libraries discussed in the previous sections (see FigfTOl-right). 
However, when using the 24 yum band, the three libraries vary 
by a factor of ~ 0.3 dex. Hence using 70 yum data provide an 
estimate of Lir nearly independent of the SED library used. 



We have seen in the previous section that the CE01 library of 
template SEDs provide a good estimate of the 70 yum emission of 
galaxies based on observed 24 yum data. Therefore, we will use 
that same library to then extrapolate from observed 70 yum to Lir. 



5. 1. Luminosity function 

The total infrared luminosity function is derived using the 
1/VWr method and the completeness correction factor derived 
from our simulations. The two GOODS fields were combined 
and then the EGS, ECDFS-O and GOODS fields were treated 
separately. Error bars were derived using Monte Carlo sim- 
ulations to take into account most of the uncertainties. For 
these Monte Carlo simulations the total infrared ^-correction 
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Fig. 9. Rest-frame 35 microns LF estimated independently in three redshift bins with the 1 /V„ mx method. The red stars are obtained 
from the GOODS-North and South fields using 70 //m sources. The dark blue and light blue stars show the LF obtained from 70 //m 
sources in the EGS and the ECDFS-O fields, respectively. The red and dark blue squares are obtained starting from 24 /jm and using 
the 24-70 yum corr elation of Figure [8] in the GOODS and EGS fields, respectively. Asterisks show the local reference taken from 
Shu pe et all ([1998) and the dotted line presents the best-fit to these data points with a double power law function (see text). In the 
intermediate redshift panel, we reproduced the best fit of the LF obtained in the two other redshift bins (i.e 0.4 < z < 0.7 dot-dash 
line and 1.0 < z < 1.3 triple-dots-dash line). The light and dark shaded area span all the solutions obtained with the^ 2 minimization 
method and compatible, within 1 <x, with the LF measured using direct 70 //m observations and the LF measured using the stacking 
analysis, respectively. The inset plot represents the evolution of the (pkme and Lk„ ee as function of redshift. Symbols and lines are as 
in Figure [6] 



uncertainties were assumed to be of 0. 1 1 dex. 

We use as a l ocal reference for the total IR LF, the 
ISanders et alj d2003l) LF derived from the analytical fit to 
the IRAS Revised Bright Galaxy Sample, i.e if> oc Zr 06 for 
log(L/L )< L knee and 4> oc Zr 2 - 2 for log(L/L G )> L knee with 
/o£io(Z.ju ee /L Q )=10.5 (z~0). 

As for the rest-frame 35 yum LF, we use the observed 24 
yum to constrain the faint end slope of the LF. For each source 
we first derived its 70 yum flux density using the 24-70 /im 
correlations that we obtained from our stacking analysis (see 
Section [4~2l . Then, as for direct 70 fim detections, we used the 
CE01 SEDs to derive Lir from the 70 /im flux density. 



We separate the values of the LF derived from observed 70 
yum (stars) and stacked 70 y um (squares) in the r esulting total IR 
LF (Figure ITTb. Using the ISanders et alj d2003l) local reference 
with fixed slopes we fit independently the total infrared LF 
derived using direct 70 yum observations (light shaded area) and 
the one derived using the stacking analysis (dark shaded area). 
We note that direct 70 //m observations allow us to sharply 
constrain the evolution of LIRGs and ULIRGs up to z ~ 1.15 
but fail to constrain the knee of the total IR LF which requires 
to rely on stacking (see light shaded area). In the regime where 
direct 70 yimi data is available, we find that the two techniques 
(starting from either 24 yum + stacking or 70 yum) provide 
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Fig. 11. Total infrared LF estimated for three re dshift bins with the \J Vmax method. Lines and symbols are the same as in Figure 
[9] Asterisks show the local reference taken from lSanders et all d2003l) and the dotted line presents the fit of the se data points with 
a dou ble power law (see text). The empty triangles and the empty circles are taken from lLe Floc'h et all d2005l) and iHuvnh et alj 
d2007l) . respectively. The light and dark shaded area spans all the parametric solutions obtained with the x 1 minimization method 
and compatible, within 1 <x, with the LF measured using direct 70 /mi observations and the LF measured using the stacking analysis 
respectively. The inset plot represents the evolution of the (/>^„ ee and Lk„ ee as function of redshift. Symbols and lines are the same as 
in Figure [6] 



consistent bright sides of the LF. The resulting parametric fits 
(using the parametric form of the local reference at all redshifts) 
are summarized in Table [5] 

The total IR LF is found to exhibit a rapid evolution with 
redshift (Fig. ITTb nearly consistent with a pure luminosity 
evolution proportional to (1 + 7) 3 - 6±0 - 4 from z ~ to z ~ 1.3. 
Marginal evidence for a decrease of 4>knee is found proportional 
to (1 + z y 0&±0 - 6 when fitting all fields, but staying consistent 
with a fiat behavior when accounting for cosmic variance (see 
Fig. HD. 

Suc h decrease of 4>knee w a s also found by previous studies 
such a s iPerez-Gonzalez et alj d2005l) for the rest-frame 12 /mi 
LF and lCaputi et alj ( 2007 ) for the rest-frame 8 //m and the total 
infrared LF. Nevertheless this decrease has to be taken with 
caution since incompleteness in the faint luminosity bins could 



mimic this evolution (Perez-Gonza lez et al .1 120051) . 

Finally, we also notice that the redshift evolution of the 
rest-frame 35 /mi and total IR LF are in agreement, i.e a pure 
luminosity evolution proportional to (1 +^) 3 - 6±0 - 5 and (1 +z) 36±0 - 4 
for the rest-frame 35 //m and total IR LF respectively. This 
consistency arises from the quasi-linearity between L35 and Ljr. 



5.2. Comparison with previous work 

As revealed in Table|6]the global evolution of the total infrared 
LF infe rred in this study is broadly consistent with previous 
studie s dLe Floc'h et al J 120051; ICaputi et af] l2007t IHuvnh et all 
|2007|) . However slight differences appear between in the LIRG 
and ULIRG regimes. In Figure [TT] we compare the LF inferred 
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Table 5. Parameter values of the 15 fan, 35 ym and total infrared LF 



Redshift 


Wavelength 




«'2 a 


Log(L fc ,„) b 


Log(<fa„ M ,) b 










Log(L ) 


Log(Mpc 




z~0 


15 fim 


-0.57 


-2.27 


9.56 ± 


0.04 


-2.73 ± 


0.07 


0.4 < z < 0.7 




-0.57 


-2.27 


10.22 i 


0.03 


-2.63 ± 


0.05 


0.7 < z < 0.7 


15 fjm 


-0.57 


-2.27 


10.57 4 


0.04 


-2.86 ± 


0.04 


1.0<z<0.7 


15 fjm 


-0.57 


-2.27 


10.79 4 


0.05 


-2.93 ± 


0.06 


z~0 


35 /im 


-0.55 


-1.95 


9.85 ± 


0.07 


-2.83 ± 


0.10 


0.4 < z < 0.7 


35 fjm 


-0.55 


-1.95 


10.73 4 


0.04 


-2.98 ± 


0.05 


0.7 < z < 1.0 


35 yum 


-0.55 


-1.95 


10.82 4 


0.04 


-2.73 ± 


0.06 


1.0 <z< 1.3 


35 yum 


-0.55 


-1.95 


11.20 4 


0.04 


-2.98 ± 


0.05 


z~0 


1R 


-0.60 


-2.20 


10.48 4 


0.02 


-2.52 ± 


0.03 


0.4 < z < 0.7 


IR 


-0.60 


-2.20 


11.194 


0.04 


-2.84 ± 


0.06 


0.7 < z < 1.0 


IR 


-0.60 


-2.20 


11.37 4 


0.03 


-2.65 ± 


0.05 


1.0<z< 1.3 


IR 


-0.60 


-2.20 


11.694 


0.06 


-2.91 ± 


0.10 



a Fixed slopes 

b These values are inferred using direct 70 fim observation and the stacking analysis 
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Fig. 10. The ^-correction between the observed 24 //m luminos- 
ity and the total infrared luminosity as a function of redshift 
(left) and the ^-correction between the observed 70 yum lumi- 
nosity and the total infrared luminosity as a function of redshift 
(right). Lines are as in Figure|5] 



in our study with the work of iLe Floc'h et a D (12005b and 
iHuvnh etalJd2007h . 

We notice that the LF inferred by ILe Floc'h etail d2005l) 
is consistent with our valu es at low lumin o sities. However, 
above log(L IR /L G ) >1 1.5, [Le Floc'h etaf] d2005l) obtained 
higher values, with a discrepancy increasing with luminosity. 
These differences, which affect the relative contributions of 
LIRGs and ULIRGs to the total infrared luminosity density, 
might be explaine d by the fact that we removed AGNs from 
our catalogs while ILe Floc'h et alj d2005l) did not. However we 
checked that keeping AGNs to derive the total infrared LF only 
reduces part of the discrepancy. Hence the main explanation of 
these differences is the method used to convert monochromatic 
lumino sities into total infrared luminosities. ILe Floc'h et alj 
(2005) used a combination of SED libraries to derive the total 
infrared luminosity from 24 fim observations while we used the 
CE01 library and either direct or stacked 70 yum observations 



0.85<z<1.15 




9 10 11 12 13 14 

Log [L, R (L e )] 



Fig. 12. Total infrared LF estimated at z ~ 1 with the l/V max 
method. Symbols are the same as in FigurefTTIex cept for the blue 
circles which represent the LF inferred by Caputi et afj d2007l) 
at z ~ 1 using 24 fim observations and the empty diamonds 
which represent the LF that we would hav e infer red using our 
24 fim observations and the Lag ache et alj d2003l) SED library 
to convert monochromatic observations into total infrared lumi- 
nosity. The shaded area represents the surface spanned by all the 
parametric solutions obtained from the x 1 minimization method 
and compatible within 1 <x. The bottom panel shows the differ- 
ences between all the observations and the best fit that we ob- 
tained with a double power law parameterization (i.e the dash 
line of the upper graphic). The solid curves are the best estimate, 
maximum and minimum values allowed by the galaxy counts at 
15,24,70,170 and 850yt/m in Le Borgne et al. (A&A, submitted). 

for the bright and faint end slopes of the LF, respectively. The 
latter should provide more robust estimates of L/« since the 
rest-frame 35 yum luminosity is produced by similar carriers as 
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Fig. 13. Evolution of the comoving number density up to z ~ 1.3 
of normal galaxies (i.e 10 7 L < L; r < 10 1 1 L ; black filled trian- 
gle), LIRGs (i.e 10 11 L < L; r < 1O 12 L ; orange filled diamond) 
and, ULIRGs (i.e 1O 12 L < L[ r ; red filled star). The green cir- 
cles represent the total number of galaxies which are above the 
detection limit of the surveys p resented here, i .e lJ' u x hm " < L[ r . 
The z ~ points are taken from lSanders et al.l d2003l) . 



the rest-fra me 24 fjm one which was proven to be a robust SFR 
indicator bv lCalzetti et all (120071) . 

We notice on Figure [TTlthat o ur LF is, within the error bars, 
in agreement with Huynh et al. (2007). Nevertheless we also 
notice that at low luminosities, their values are slightly lower 
than ours. We interpret this discrepancy as a difference in the 
treatment of incompleteness. Our multi-wavelengths catalogs 
being computed with a new extraction technique based on prior 
positions, we are able to extract blended sources and hence 
reduce incompleteness corrections. We conclude that at low 
luminosities the LFs measure d in this study sho uld be more 
robust than the ones derived by Huvnh etalJ (120071) . 

The total IR LF derived here is noticeably differ ent 
from the one previously obtained by ICaputi et al.l d2007l) at 
Z ~ 1 (see Figure [T2l who used 24 fim images in the CDFS. 
There are two possible explanations for this difference: it could 
either result from cosmic variance or from the method that they 
used to derive L/« and which does not take into account 70 fim 
observations as in our study. To exclude the first possibility, 
we computed the total IR LF from our EGS 24 yum catalog. 
Depending on the SED library that we used, we either find an 
IR LF overlapping ours (using the CE01 SEDs) or matching 
the ICaputi et al.l (120071) one (using the vLy 4/ ™-LrR correlation 
of Bavouzet et al. 2008 as in Caputi et al. 2007). This shows 
that the difference between the two studies does not come from 
cosmic variance but rather from the choice of the SED library. 
As we showed earlier, the CE01 SED turned out to fit better the 
24 versus 70 fim correlation down to the faintest luminosities 
using stacking. Since the 70 fim band provides a more accurate 
estimate of the total infrared luminosity than the 24 /im one and 
since the LDP library does not reproduce the observed 24 - 70 
jjm correlation at z — 1 (see Figure |8j, we conclude that our 
inferred LFs should be more robust. 



6. Discussion 

The redshift evolution of the comoving number density of 
LIRGs and ULIRGs is derived by integrating the infrared LF in 
each redshift bin, 0.4 < z < 0.7, 0.7 < z < 1.0 and 1.0 < z < 1.3. 
The number density of LIRGs and ULIRGs appear to be 40 and 
100 times larger at z ~ 1 than in the local Universe, respectively 
(see Figure ITS!. Hence we confirm the rapid evolution of these 
extremely luminous dusty galaxies over the last 8 billion years. 
The evolution is at a slightly more m ode rate pace than previous 
studies since iLe Floc'h et~atl (120051) and ICharv & Elbazl (1200 lb 
found a factor 70 for LIRGs. This correction towards a larger 
contribution of the remaining lower luminosity galaxies, that 
we name here "normal" galaxies, comes from the faint part 
of the LF that we were able to constrain through our stacking 
analysis. It must be noted that at the highest redshifts, we are 
only sensitive to the LIRG and ULIRG population Thus the 
extrapolation to fainter luminosities (illustrated by the green 
filled circles in Fig.fTJt assumes the same faint end slope for the 
LF seen at lower redshifts and could be uncertain. Finally, we 
note that the number density of ULIRGs at z ~ 1.15 is nearly 
equal to that of LIRGs at z ~ 0. In order to detect 10 ULIRGs at 
z ~ 1+0.2 (892 Mpc 3 /arcmin 2 ), one would need to cover more 
than 3 GOODS fields. 

In Figure [14] we compare the comoving IR luminosity 
density (or equivalently SFR density due to obscured star 
formation, using Eq. |4]i produce d by normal galaxies, LIRGs 
and ULIRGs to that obtained by ILe Floc'h etaf] d2005l) . Note 

that X-ray AGN have already been removed fro m the LF. 

W e also notice that at z ~ values quoted in lLe Floc'h et al.l 
(120051) need to be corrected by a factor 2.5 (empty triangles of 
Figure [T4l. Th e total infrared LF tak en as their local reference 
was derived in Sa nders et al.l d2003l) by bins of magnitude but 
used in ILe Floc'h et al.l d2005l) as if it was derived by bins of 
luminosity. This error does n o t chan ge the LFs derived at higher 
redshift by ILe Floc'h et al.1 d2005l) but could have affected 
their estimates of ar<j and and hence their estimates of the 
evolution of the infrared luminosity density. The discrepancy 
fo und between the ad der ived in this study and the one derived 
in ILe Floc'h et al. ( 2005J) can be partly due to this error since 
starting from an erroneous local LF would yield to an overesti- 
mate of a,i. 

As observed in Figure[14]we confirm the transition between 
a regime dominated by normal galaxies at low redshift and 
by LIRGs at high redshift. The transition occur r ed at z ~ 0.9 
compared to z ~ 0.7 derived bv lLe Floc'h et al] d2005l) . While 
LIRGs and ULIRGs altogether produce less than 2 % of the 
present-day bolometric luminosity density of galaxies, they 
produced half of it 8 billion years ago. Indeed, the contribution 
of unobscured UV light to the comoving SFR density was 
lower than that of LIRGs/ULIRGs by a factor of about ~3 (see 
Fi g . |T3T> ■ Although the luminosity density of ULIRGs was much 
larger at z ~ 1 than today, their contribution to the total IR 
energy density was still rathe r low ( i.e ~ 10%). We find it to 
be lower than ILe Floc'h et al.l d2005l) by a fact or of ~ 2 which 
may p artly be explained by the fact that unlike ILe Floc'h et al.l 
(120051) we have removed the contribution of X-ray AGN which 
becomes large among ULIRGs. Furthermore, by using the 
observed 70 /im data instead of extrapolating from 24 fim, we 
have reduced the uncertainty in the estimates of the bolometric 
luminosity. 

In order to get a complete census on the SFR history we 
need to take into account the contribution of unobscured UV 
light. The unobscured SFR density (dot line in Figure |T3T > is 
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Table 6. Parameters leading the total infrared LF evolution 



Authors 


Functional Form" 


Wavelength 






a? 


This study 

Le Floc'h et al. (2005) 
Caouti et al. (2007) 

Huvnh et al. (2007) 
Hopkins (2004) 


DPL 
DE 
DE 

DPL 


24 & 70 fjm 
24 fim 
24 fim 

70 fim 
SFR compilation 


z< 1 
z> 1 


3.6 ± 0.4 

t 9+0.7 
•'-0.2 

3.5 ± 0.4 
2.2 + 0.5 

9 -7O+0.34 
L /8 -0.32 

2.7 + 0.6 


-0.8 + 0.6 

7 +0 ' 2 

-0.6 

-0.7 + 0.1 
-3.9+ 1.0 

0.0 
0.15 + 0.6 



DPL: Double Power Law; DE: Double Exponential 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Redshift 



Fig. 14. Evolution of the comoving IR energy density up to z ~ 1.3 (striped area) and the relative contribution of normal galaxies 
(i.e 7 < log(Lj r ) < 11; yellow filled area), LIRGs (i.e 11 < log(Lj r ) < 12; orange filled area) and, ULIRGs (i.e 12 < log(Li r ); red 
filled area). The areas w ere defined using all th e solutions compatible within 1 cr with the total infrared LF. Black triangles represent 
the results obtained by lLe Floc'h et alJ (120051) for the global evolution of the comoving energy density (solid line) and the relative 
contribution of normal galax ies (triple dot dash li ne), LIRGs (dashed line) and, ULIRGs (dot dash line). At z ~ empty black 
triangles represent results of lLe Floc'h et alJ d2005l) corrected by a factor 2.5 (see text). 



taken from a GALEX study ( Schimi novich et ail 120051) which 
yield an evolution following a (1+z) 2 5 law and a local star 
formation density of psfr(z ~0)=4 x 1(T 3 M yr _1 Mpc~ 3 . The 
total SFR density (dashed line in Figure |T3T> was then defined as 
the sum of the unobscured SFR density traced by the direct UV 
light and the dusty SFR density traced by the infrared emission. 
We find that the relative contribution of unobscured UV light 
to the cosmic SFR density evolves nearly in parallel with the 
total one and accounts for ~2Q % of the local SFR density and 
-20% of the p SFR (z ~1.15)=0.15 Moyr-'Mpc" 3 . Globally, the 
cosmic star-formation history that we derived using data at 24 
and 70 fim is consistent with the combination of indicators, 
eit her unobscured or co r recte d for dust extinction, as compiled 
by [Hopkins & Beacoml d2006l) . We also notice a very good 
agreement between the cosmic star -formation history de rived 
in our work and the ones derived by Seym our et alJ d2008l) and 
Smo lcic et alJ (2009) using deep radio observations. Especially, 



it is very encouraging to see that the star formation history 
for ULIRGs obtained using deep radio data is in very good 
agreement with our independent measurements. 



7. Conclusion 

We make use of the deep 24 fim and 70 fim surveys in the Great 
Observatories Origins Deep Survey (GOODS) and Extended 
Groth Strip (EGS) fields to characterize the evolution of the 
comoving star-formation rate density at 0.4 < z < 1.3 which 
corresponds to 2/3 of cosmic time. The data span a total area of 
~800 arcmin 2 and a comoving volume of 1.4x1 6 Mpc 3 . Thanks 
to the wealth of multiwavelength data, ~80% of the sources are 
associated with either a photometric or spectroscopic redshift. 

The exquisite depth of the data results in a high surface 
density of sources. In order to alleviate the effects of confusion, 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Redshift 



Fig. 15. Evolution of the comoving star formation rate density up to z ~ 1.3 (striped area) and the relative contribution of normal 
galaxies (i.e 7 < log(Li r ) < 11; yellow filled area), LIRGs (i.e 1 1 < log(L; r ) < 12; orange filled area) and, ULIRGs (i.e 12 < log(Li r ); 
red filled area). The areas were defined using all the solutions compatible with the total infrared LF within 1 cr. The solid line 
represents the best fit to the total SFR density. The dottes line represents the SFR measured using the UV light not corrected from dust 
extinction. The dashed line represents the total SFR density defined as the sum of the SFR density estimated b y the infrared and of 
the SF R density obtained from the UV light uncorrected of dust extinction. Light blue diamonds are taken from Hop kins & Beacon] 
(200 6j) and represent the SFR densities estimates using various estimators. Dark blue triangles represent the SFR de nsity estimated 
by ISeymour etaD (2008) using deep radio observations. Dark blue circles represent the SFR density estimated by ISmolcic et alj 
mm using deep 20 cm observations and dark blue squares represent the relative contribution of ULIRGs to this SFR density. 



we have developed a technique which uses the prior positional 
information of sources at shorter wavelengths to extract the flux 
density of sources at 24 fim and 70 fim. We have undertaken 
extensive simulations to quantify flux errors and biases which 
might be associated with this technique. The relative depth of 
the 24 and 70 fim survey results in only 7% of 24 fim sources 
being detected at 70 fim. 

By associating the redshifts with the sources in the mid- 
and far-infrared catalogs, we are able to measure the rest-frame 
15 fim and 35 fim luminosities of galaxies out to z ~ 1.3. We use 
a stacking procedure to measure the evolution of the 24 to 70 fim 
luminosity ratio for sources which are undetected at 70 /urn. We 
find that the average observed 24yum and 70 fim properties of 
both the detected and stacked sources are in excellent agreement 
with correlations derived from local galaxy samples. However, 
the 24/70 micron flux ratios of individual galaxies do show a 
significant scatter of 0.2 dex. 

By complementing the direct detections with the stacking 
analysis, we are able to measure a 15 fim, 35 fim and total 
infrared luminosity function for star forming galaxies between 
0.4 < z < 1.3. We measure both the bright end and faint 
end slope of the luminosity functions, which span a factor of 
-100 in luminosity (10 n L < L IR < 10 13 L ), a substantial 
improvement over previous works. We find no evidence for a 
change in slope of the double power-law used to characterize 
the local IRLF of galaxies. However, we confirm the strong 



evolution in the number density of infrared luminous galaxies 
between < z < 1 . 

LIRGs and ULIRGs evolve at about the same rate over 
this redshift, having increased by a factor of ~40 and ~100 in 
number density, respectively. LIRGs contribute about ~50% to 
the comoving star-formation rate density at z ~ 1, comparable 
to the contribution of lower luminosity star-forming galaxies. 
The infrared luminosity density is a factor of five higher than 
the UV luminosity density at z ~ 1 confirming that the bulk of 
the star-formation is dust obscured. 
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